64 research outputs found
Collisional dark matter density profiles around supermassive black holes
We solve the spherically symmetric time dependent relativistic Euler
equations on a Schwarzschild background space-time for a perfect fluid, where
the perfect fluid models the dark matter and the space-time background is that
of a non-rotating supermassive black hole. We consider the fluid obeys an ideal
gas equation of state as a simple model of dark matter with pressure. Assuming
out of equilibrium initial conditions we search for late-time attractor type of
solutions, which we found to show a constant accretion rate for the non-zero
pressure case, that is, the pressure itself suffices to produce stationary
accretion regimes. We then analyze the resulting density profile of such
late-time solutions with the function . For different values of
the adiabatic index we find different slopes of the density profile, and we
study such profile in two regions: a region one near the black hole, located
from the horizon up to 50 and a region two from up to , which for a black hole of corresponds to pc. The profile depends on the adiabatic index or equivalently on the
pressure of the fluid and our findings are as follows: in the near region the
density profile shows values of and in the limit of the
pressure-less case ; on the other hand, in region two,
the value of in all the cases we studied. If these results are to
be applied to the dark matter problem, the conclusion is that, in the limit of
pressure-less gas the density profile is cuspy only near the black hole and
approaches a non-cuspy profile at bigger scales within 1pc. These results show
on the one hand that pressure suffices to provide flat density profiles of dark
matter and on the other hand show that the presence of a central black hole
does not distort the density profile of dark matter at scales of 0.1pc.Comment: 7 pages, 8 eps figures, accepted for publication in MNRA
Relativistic Hydrodynamics around Black Holes and Horizon Adapted Coordinate Systems
Despite the fact that the Schwarzschild and Kerr solutions for the Einstein
equations, when written in standard Schwarzschild and Boyer-Lindquist
coordinates, present coordinate singularities, all numerical studies of
accretion flows onto collapsed objects have been widely using them over the
years. This approach introduces conceptual and practical complications in
places where a smooth solution should be guaranteed, i.e., at the gravitational
radius. In the present paper, we propose an alternative way of solving the
general relativistic hydrodynamic equations in background (fixed) black hole
spacetimes. We identify classes of coordinates in which the (possibly rotating)
black hole metric is free of coordinate singularities at the horizon,
independent of time, and admits a spacelike decomposition. In the spherically
symmetric, non-rotating case, we re-derive exact solutions for dust and perfect
fluid accretion in Eddington-Finkelstein coordinates, and compare with
numerical hydrodynamic integrations. We perform representative axisymmetric
computations. These demonstrations suggest that the use of those coordinate
systems carries significant improvements over the standard approach, especially
for higher dimensional studies.Comment: 10 pages, 4 postscript figures, accepted for publication in Phys.
Rev.
An improved formulation of the relativistic hydrodynamics equations in 2D Cartesian coordinates
A number of astrophysical scenarios possess and preserve an overall
cylindrical symmetry also when undergoing a catastrophic and nonlinear
evolution. Exploiting such a symmetry, these processes can be studied through
numerical-relativity simulations at smaller computational costs and at
considerably larger spatial resolutions. We here present a new
flux-conservative formulation of the relativistic hydrodynamics equations in
cylindrical coordinates. By rearranging those terms in the equations which are
the sources of the largest numerical errors, the new formulation yields a
global truncation error which is one or more orders of magnitude smaller than
those of alternative and commonly used formulations. We illustrate this through
a series of numerical tests involving the evolution of oscillating spherical
and rotating stars, as well as shock-tube tests.Comment: 19 pages, 9 figure
Accurate evolutions of inspiralling neutron-star binaries: assessment of the truncation error
We have recently presented an investigation in full general relativity of the
dynamics and gravitational-wave emission from binary neutron stars which
inspiral and merge, producing a black hole surrounded by a torus (see
arXiv:0804.0594). We here discuss in more detail the convergence properties of
the results presented in arXiv:0804.0594 and, in particular, the deterioration
of the convergence rate at the merger and during the survival of the merged
object, when strong shocks are formed and turbulence develops. We also show
that physically reasonable and numerically convergent results obtained at
low-resolution suffer however from large truncation errors and hence are of
little physical use. We summarize our findings in an "error budget", which
includes the different sources of possible inaccuracies we have investigated
and provides a first quantitative assessment of the precision in the modelling
of compact fluid binaries.Comment: 13 pages, 5 figures. Minor changes to match published version. Added
figure 5 right pane
On the Shear Instability in Relativistic Neutron Stars
We present new results on instabilities in rapidly and differentially
rotating neutron stars. We model the stars in full general relativity and
describe the stellar matter adopting a cold realistic equation of state based
on the unified SLy prescription. We provide evidence that rapidly and
differentially rotating stars that are below the expected threshold for the
dynamical bar-mode instability, beta_c = T/|W| ~ 0.25, do nevertheless develop
a shear instability on a dynamical timescale and for a wide range of values of
beta. This class of instability, which has so far been found only for small
values of beta and with very small growth rates, is therefore more generic than
previously found and potentially more effective in producing strong sources of
gravitational waves. Overall, our findings support the phenomenological
predictions made by Watts, Andersson and Jones on the nature of the low-T/|W|.Comment: 20 pages; accepted to the Classical and Quantum Gravity special issue
for MICRA200
Gravitational waves from axisymmetrically oscillating neutron stars in general relativistic simulations
Gravitational waves from oscillating neutron stars in axial symmetry are
studied performing numerical simulations in full general relativity. Neutron
stars are modeled by a polytropic equation of state for simplicity. A
gauge-invariant wave extraction method as well as a quadrupole formula are
adopted for computation of gravitational waves. It is found that the
gauge-invariant variables systematically contain numerical errors generated
near the outer boundaries in the present axisymmetric computation. We clarify
their origin, and illustrate it possible to eliminate the dominant part of the
systematic errors. The best corrected waveforms for oscillating and rotating
stars currently contain errors of magnitude in the local wave
zone. Comparing the waveforms obtained by the gauge-invariant technique with
those by the quadrupole formula, it is shown that the quadrupole formula yields
approximate gravitational waveforms besides a systematic underestimation of the
amplitude of where and denote the mass and the radius of
neutron stars. However, the wave phase and modulation of the amplitude can be
computed accurately. This indicates that the quadrupole formula is a useful
tool for studying gravitational waves from rotating stellar core collapse to a
neutron star in fully general relativistic simulations. Properties of the
gravitational waveforms from the oscillating and rigidly rotating neutron stars
are also addressed paying attention to the oscillation associated with
fundamental modes
Relativistic hydrodynamics on spacelike and null surfaces: Formalism and computations of spherically symmetric spacetimes
We introduce a formulation of Eulerian general relativistic hydrodynamics
which is applicable for (perfect) fluid data prescribed on either spacelike or
null hypersurfaces. Simple explicit expressions for the characteristic speeds
and fields are derived in the general case. A complete implementation of the
formalism is developed in the case of spherical symmetry. The algorithm is
tested in a number of different situations, predisposing for a range of
possible applications. We consider the Riemann problem for a polytropic gas,
with initial data given on a retarded/advanced time slice of Minkowski
spacetime. We compute perfect fluid accretion onto a Schwarzschild black hole
spacetime using ingoing null Eddington-Finkelstein coordinates. Tests of fluid
evolution on dynamic background include constant density and TOV stars sliced
along the radial null cones. Finally, we consider the accretion of
self-gravitating matter onto a central black hole and the ensuing increase in
the mass of the black hole horizon.Comment: 23 pages, 13 figures, submitted to Phys. Rev.
Accurate evolutions of unequal-mass neutron-star binaries: properties of the torus and short GRB engines
We present new results from accurate and fully general-relativistic
simulations of the coalescence of unmagnetized binary neutron stars with
various mass ratios. The evolution of the stars is followed through the
inspiral phase, the merger and prompt collapse to a black hole, up until the
appearance of a thick accretion disk, which is studied as it enters and remains
in a regime of quasi-steady accretion. Although a simple ideal-fluid equation
of state with \Gamma=2 is used, this work presents a systematic study within a
fully general relativistic framework of the properties of the resulting
black-hole--torus system produced by the merger of unequal-mass binaries. More
specifically, we show that: (1) The mass of the torus increases considerably
with the mass asymmetry and equal-mass binaries do not produce significant tori
if they have a total baryonic mass M_tot >~ 3.7 M_sun; (2) Tori with masses
M_tor ~ 0.2 M_sun are measured for binaries with M_tot ~ 3.4 M_sun and mass
ratios q ~ 0.75-0.85; (3) The mass of the torus can be estimated by the simple
expression M_tor(q, M_tot) = [c_1 (1-q) + c_2](M_max-M_tot), involving the
maximum mass for the binaries and coefficients constrained from the
simulations, and suggesting that the tori can have masses as large as M_tor ~
0.35 M_sun for M_tot ~ 2.8 M_sun and q ~ 0.75-0.85; (4) Using a novel technique
to analyze the evolution of the tori we find no evidence for the onset of
non-axisymmetric instabilities and that very little, if any, of their mass is
unbound; (5) Finally, for all the binaries considered we compute the complete
gravitational waveforms and the recoils imparted to the black holes, discussing
the prospects of detection of these sources for a number of present and future
detectors.Comment: 35 pages; small changes to match the published versio
Axisymmetric general relativistic hydrodynamics: Long-term evolution of neutron stars and stellar collapse to neutron stars and black holes
We report a new implementation for axisymmetric simulation in full general
relativity. In this implementation, the Einstein equations are solved using the
Nakamura-Shibata formulation with the so-called cartoon method to impose an
axisymmetric boundary condition, and the general relativistic hydrodynamic
equations are solved using a high-resolution shock-capturing scheme based on an
approximate Riemann solver. As tests, we performed the following simulations:
(i) long-term evolution of non-rotating and rapidly rotating neutron stars,
(ii) long-term evolution of neutron stars of a high-amplitude damping
oscillation accompanied with shock formation, (iii) collapse of unstable
neutron stars to black holes, and (iv) stellar collapses to neutron stars. The
tests (i)--(iii) were carried out with the -law equation of state, and
the test (iv) with a more realistic parametric equation of state for
high-density matter. We found that this new implementation works very well: It
is possible to perform the simulations for stable neutron stars for more than
10 dynamical time scales, to capture strong shocks formed at stellar core
collapses, and to accurately compute the mass of black holes formed after the
collapse and subsequent accretion. In conclusion, this implementation is robust
enough to apply to astrophysical problems such as stellar core collapse of
massive stars to a neutron star and black hole, phase transition of a neutron
star to a high-density star, and accretion-induced collapse of a neutron star
to a black hole. The result for the first simulation of stellar core collapse
to a neutron star started from a realistic initial condition is also presented.Comment: 28 pages, to appear in PRD 67, 0440XX (2003
A New Open-Source Code for Spherically-Symmetric Stellar Collapse to Neutron Stars and Black Holes
We present the new open-source spherically-symmetric general-relativistic
(GR) hydrodynamics code GR1D. It is based on the Eulerian formulation of GR
hydrodynamics (GRHD) put forth by Romero-Ibanez-Gourgoulhon and employs
radial-gauge, polar-slicing coordinates in which the 3+1 equations simplify
substantially. We discretize the GRHD equations with a finite-volume scheme,
employing piecewise-parabolic reconstruction and an approximate Riemann solver.
GR1D is intended for the simulation of stellar collapse to neutron stars and
black holes and will also serve as a testbed for modeling technology to be
incorporated in multi-D GR codes. Its GRHD part is coupled to various
finite-temperature microphysical equations of state in tabulated form that we
make available with GR1D. An approximate deleptonization scheme for the
collapse phase and a neutrino-leakage/heating scheme for the postbounce epoch
are included and described. We also derive the equations for effective rotation
in 1D and implement them in GR1D. We present an array of standard test
calculations and also show how simple analytic equations of state in
combination with presupernova models from stellar evolutionary calculations can
be used to study qualitative aspects of black hole formation in failing
rotating core-collapse supernovae. In addition, we present a simulation with
microphysical EOS and neutrino leakage/heating of a failing core-collapse
supernova and black hole formation in a presupernova model of a 40 solar mass
zero-age main-sequence star. We find good agreement on the time of black hole
formation (within 20%) and last stable protoneutron star mass (within 10%) with
predictions from simulations with full Boltzmann neutrino radiation
hydrodynamics.Comment: 25 pages, 6 figures, 2 appendices. Accepted for publication to the
Classical and Quantum Gravity special issue for MICRA2009. Code may be
downloaded from http://www.stellarcollapse.org Update: corrected title, small
modifications suggested by the referees, added source term derivation in
appendix
- …